Gigantic Enhancement of Magnet o-Chiral Effect in Photonic Crystals 
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We theoretically propose a method to enhance dramatically a magneto-chiral(MC) effect by using 
the photonic crystals composed of a multiferroic material. The MC effect, the directional birefrin- 
gence even for unpolarized light, is so small that it has been difficult to observe experimentally. 
Two kinds of periodic structures are investigated; (a) a multilayer and (b) a stripe composed of 
a magneto-chiral material and air. In both cases, the difference in reflectivity between different 
magnetization directions is enhanced by a factor of hundreds compared with a bulk material. 
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Optical effects in media with broken symmetry have 
been of great interest from both scientific and ap- 
plicational viewpoints. Broken inversion symmetry X 
and time-reversal symmetry T give rise to polarization- 
dependent properties; natural optical activity and 
magneto-optical effects, respectively. Polarization- 
independent properties require breaking of both X and 
T . When both X and T are simultaneously broken, we 
have directional birefringence even for unpolarized light, 
which is non-reciprocal. Such directional birefringence 
is conventionally called a magneto-chiral (MC) effect 1 
or an optical magnetoelectric(OME) effect 2 , if the mi- 
croscopic structure of the medium is helical or polar, re- 
spectively. Nevertheless, as the polarization-independent 
non-reciprocal effects, the MC effect and OME effect are 
quite similar, and we call both the effects "MC effects" 
for simplicity. The MC effect has been mostly studied in 
chiral media in external fields 1 ' 3 . On the other hand, only 
recently has it been investigated in condensed materials 
without external fields. Such materials should be multi- 
ferroic, i.e., without X and T. Multiferroic materials such 
as GaFe03 exhibit the spontaneous MC effect even in the 
absence of the external fields 2 ' 4 ' 5 . Optical properties in 
multiferroics can be described by a toroidal moment de- 
fined by T = J2i TiXSi ~ Px M, where Vi is the displace- 
ment from the center position of atoms, S{ is a magnetic 
moment at i-th site, P is the spontaneous electric polar- 
ization, and M is the magnetization, respectively 6 . The 
dielectric function depends on whether light propagation 
is parallel or antiparallel to the toroidal moment. 

Unfortunately the MC effect is usually too small to be 
observed experimentally. In this letter we theoretically 
propose a method to enhance the MC effect by using a 
photonic crystal. A photonic crystal is a periodic array 
of dielectric materials 7 ' 8 ' 9 . Light propagation in a pho- 
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tonic crystal described by a wave equation is analogous 
to electron motion in a solid. In this analogy the wave 
equation with a spatially modulated dielectric constant 
corresponds to the Schrodinger equation with a poten- 
tial. However, most studies on photonic crystals do not 
utilize this analogy to electronic systems; instead they 
aim only to realize the localization of light or high de- 
values. In condensed materials, magnetism has been one 
of the central issues, while much less attention has been 
paid to magnetic photonic crystals 10 . In light of rich 
physics and applications arising from magnetism in con- 
densed materials, magnetic photonic crystals should have 
potential importance both for fundamental physics and 
for application. Thus it is highly desirable to construct a 
theory describing light propagation in materials without 
time-reversal symmetry. 

We theoretically found that the difference in reflectiv- 
ity for opposite directions of the toroidal moment is mag- 
nified by hundreds of times in photonic crystals than that 
in a bulk. The enhancement is due to the Bragg reflec- 
tion and a photonic band gap, which are absent in a bulk 
or a Fabry-Perot cavity. In this respect photonic crystals 
are of great promise for better interferometers in place of 
Fabry-Perot cavities. 

We solve the following wave equation derived from 
Maxwell equations, x [V x E(r)] = (^) 2 £(f), 

where e(f) is a dielectric function and E is an electric 
field vector. The direction of the toroidal moment is 
postulated to be parallel to the x axis, and the MC 
effect can be introduced in the dielectric function of 
the medium 1 as S\ = e\ + ak x: where 6® is the static 
dielectric constant, a is a strength of the MC effect 
and k x is an operator representing the x component of 
the unit wavevector. We set the dielectric constants of 
the two media as e\ = 15, a = 10 -6 and £2 = l(air), 
and consider the normal incidence throughout this 
letter. In order to solve the equation above, we use a 
transfer matrix method for multilayer structures shown 
in Fig. 1(a), and a plane- wave expansion method 11 for 
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stripe structures shown in Fig. (b). The polarization 
change, i.e., magneto-optical effects and natural optical 
activity, are neglected throughout this letter, and we 
only focus on the MC effect, which is polarization- 
independent. We calculate frequency dependence of 
the difference in reflectivity for opposite directions 
of the toroidal moments. It should be noted that 
the difference is zero in a conventional material, and 
non-zero only in a magneto-chiral medium. First 
of all, we calculate reflectivity without the periodic 
modulation of e(f). Let us consider light reflection of 
a magneto-chiral material with the normal incidence; 
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where the indices j and j denote whether f is parallel 
or antiparallel to the x-axis. This value should be 
compared with that in a periodic structure. A MC 
effect in a material is originated from the spin-orbit 
interaction that is nothing but a relativistic effect, 
and is difficult to control microscopically. In con- 
trast, we can enhance the effect by the artificial periodic 
modulation with the size of the order of light wavelength. 

(a) Magneto-chiral multilayers 

The multilayered structure (16 layers) composed of the 
MC medium and air is shown in Fig. 1(a). The ratio of 
the thickness between the MC medium and air is put 
to be 0.2 : 0.8. We solve the wave equation by means 
of the transfer matrix method. The boundary conditions 
and the phase change through the propagation can be de- 
scribed as matrices, and we can solve the equation simply 
by multiplying the matrices. 

We calculate reflectivity R^^ where the arrows rep- 
resent the directions of T whether it is parallel or an- 
tiparallel to the x axis. Figure 2 shows the frequency 
dependence of the difference in reflectivities for opposite 
directions of T. Here we introduced small absorption. 
There are two distinct behaviors in the figure. One is the 
region characterized by "gap" in the figure corresponds 
to the photonic band gap. The other is the region where 
Rl —R^ oscillates, corresponding to a Fabry-Perot cavity 
that has nothing to do with the periodic structure. At 
the edge of the gap, we have a remarkable enhancement. 
This enhancement is interpreted as follows. The different 
toroidal moment direction has the different refractive in- 
dex, leading to different reflectivity. Hence the band gaps 
for the opposite toroidal moment directions are at differ- 
ent frequency regions. At the gap edge, it can happen 
that the frequency is in the band gap for T \\ e x , and in 
the conduction band for T \\ —e Xl because we have the 
nonreciprocal dispersion uo{k) ^ uj{—k). Therefore the 
difference in reflectivities between the two toroidal mo- 
ment directions is significantly magnified. In other words, 
this difference in reflectivities behaves qualitatively sim- 
ilar to t^j, and the differential is singular at the edge 
of the gap, resulting in the enhancement. The magnifi- 
cation at the gap edge is robust under small changes of 



incident angles and toroidal moment direction, whereas 
the values in the Fabry-Perot region are fragile against 
these changes. We note that if there is no absorption in 
a MC medium, the value of R} — R^ at the gap edge is 
the order of 10 -5 , an order of magnitude larger. 



(b) Magneto-chiral stripes 

In the stripe structures shown in Fig. 1(b), reflectivity 
changes as the direction of the toroidal moment is re- 
versed. We calculate R[ — R[, which is the difference 
between reflectivities of the first-order Bragg reflected 
waves, denoted by the subscript 1. We solve the equation 
by Fourier transformation in the region < z < L, 

oo 

e= J2 J2 A ™ e * knXe ^ z > t 1 ) 

n=0,±l rn=0 

where k n = ^ 21 , L is the thickness of the 
stripe, and we have a sinusoidal modulation 
i£y = Wi 1 + ^) + Wi 1 - ^cos^, where 
a a periodicity along the x direction. It should be 
noted that the structure shown in Fig. 1(b) is only a 
schematic one. In the calculation, we use the sinusoidal 
dielectric modulation instead of the rectangular one, 
because the higher order Fourier components give no 
essential contribution to the results. We choose the 
dielectric constant of the substrate as £3 = 15. We 
take 606 Fourier components, n = —1,0,1, m = 0-200, 
which are enough for convergence. Figure 3 represents 
the frequency dependence of the difference in intensities 
of Bragg reflected waves for L = 0.5a, a, and 2a. In 
the region ^ < 1.0, the wavenumber is imaginary, 
implying that the reflected light is an evanescent wave. 
Therefore we should focus on the region ^ > 1.0. In 
the region just above ^ = 1.0, the reflection angle is 
close to 90°, and it is experimentally difficult to detect 
the beam. Let us consider the case L = a at first. We 
find an enhancement at ^ ~ 1.3 corresponding to the 
Bragg reflection angle 9 = 63°, which can be observed 
experimentally. Moreover, even away from ^ c± 1.3, 

R[ — R[ is of the order of 10 -6 , which is two orders 
of magnitude larger than that in a bulk. Although 
there is no band gap in this system, the optical effect 
is remarkably enhanced when the wavelength of light is 
comparable to the periodicity of the photonic crystals. 
When the thickness of the stripe changes, L = 0.5a and 
2a, the order of rI — R[ remains the same. It is because 
this enhancement is due to the periodicity along the x 
direction and has nothing to do with the value of the 
thickness. 

Let us discuss the relation and the difference between 
(a) and (b) in Fig. 1. In the MC multilayers, a peri- 
odic structure and the non-reciprocal dispersion uo(k) ^ 
uo(—k) play essential roles. On the other hand, in the MC 
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FIG. 1: Schematic illustration of (a) a magneto- chiral multi- 
layer and (b) a stripe structure whose constituent materials 
are multiferroics and air. The vector G in (b) denotes a re- 
ciprocal lattice vector. 

stripes without any band gap, a periodic structure brings 
about multiple scatterings, and only the waves that sat- 
isfy the Bragg conditions can be diffracted. In both cases, 
w — is roughly a differential of R with respect to uj, 
where R} is or R[. The multiple scatterings induce 
the oscillation of R and the large resulting in a re- 
markable enhancement of w — R^. In order to magnify 
this effect, it is necessary to array the materials paral- 
lel to the toroidal moment direction. It should be noted 
that the enhancement factor of hundreds is not specific 
to these cases (a = 10 -6 ) but valid for another value of 
a. 

In conclusion, we have investigated an enhancement 
of the magneto-chiral effect in photonic crystals. The 
difference in intensities of two waves for different toroidal 
moment directions is hundreds of times enhanced from 
that in a bulk. The direction of the toroidal moment can 
be reversed by changing the directions of the polarization 
or the magnetization by external fields. Our theory is 
in principle applicable for a wide range of materials, by 
changing dielectric constants and thickness of the layers 
at one's disposal. 
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FIG. 2: Frequency dependence of R} — for the normal 
incidence. The small absorption is introduced as e\ — 15 + 
O.li. 
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FIG. 3: Frequency dependence of R\ — R[ for the normal 
incidence. 
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